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ABSTRACT 

It is argued that there is a linear correlation between star formation rate 
(SFR) and accretion rate for normal bright active galactic nuclei (AGNs). How- 
ever, it is still unclear whether this correlation holds for LINERs, of which the 
accretion rates are relatively lower than those of normal bright AGNs. The ra- 
diatively inefficient accretion flows (RIAFs) are believed to be present in these 
LINERs. In this work, we derive accretion rates for a sample of LINERs from 
their hard X-ray luminosities based on spectral calculations for RIAFs. We find 
that LINERs follow the same correlation between star formation rate and ac- 
cretion rate defined by normal bright AGNs, when reasonable parameters are 
adopted for RIAFs. It means that the gases feed the black hole and star forma- 
tion in these low-luminosity LINERs may follow the same way as that in normal 
bright AGNs, which is roughly consistent with recent numerical simulations on 
quasar evolution. 

Subject headings: galaxies: active — galaxies: nuclei — X-rays: galaxies — accre- 
tion, accretion disks 



1. Introduction 

Black hole (BH) accretion is thought to power AGNs, and the UV/optical bump ob- 
served in bright quasars is naturally interpreted as blackbody emission from standard thin 
accretion disks (e.g., Sun, & Malkan 1989). The AGN activity may be switched off when the 
gas near the black hole is exhausted (see Narayan 2002, for a recent review and references 
therein). Most nearby galactic nuclei are much less active and show very different properties 
with bright quasars, such as low ionization, lack "UV bump" etc (see Ho 2005, for a review). 
The accretion mode in these nearby galaxies may be different from that in powerful AGNs. 
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Narayan & Yi (1994) proposed that the standard accretion disk should transit to a RIAF 
when the accretion rate m(=M /Mftdd) declines below a critical value m cr i t within a certain 
transition radius (e.g., Esin, McClintock, & Narayan 1997; Yuan & Narayan 2004; Lu et 
al. 2004). RIAFs are optically thin, geometrically thick and very hot, which are supposed 
to be present in many low-luminosity AGNs (e.g., Lasota et al. 1996; Quataert et al. 1999; 
Gammie et al. 1999; Yuan, et al. 2002; Yuan & Narayan 2004; Nemmen et al. 2005) and our 
Galactic center Sgr A* (e.g., Narayan & Yi 1995). 

A very low nuclear-luminosity class of low- ionization nuclear emission-line region (LINER) 
galaxies were identified by Heckman (1980), in which are found approximately 2/3 low lu- 
minosity AGNs (Ho, Filippenko, & Sargent 1997a). An important fraction of LINERs are 
clearly weak manifestations of quasar-like phenomena, as demonstrated by the presence of 
broad Ha emission, which are almost certainly accretion-powered (Ho et al. 1997b). These 
LINERs can be identified with quiescent BH remnants from the quasar era. In the present 
epoch, the supply of the gas available for powering central engines is much curtailed. The 
low mass accretion rates inferred for many LINERs may suggest that these sources represent 
the "missing link" between powerful AGNs and normal galaxies as our own. 

The recent discovery of the correlation between BH mass and stellar velocity dispersion 
in nearby galaxies (e.g., Gebhardt et al. 2000; Ferrarese & Merritt 2000) demonstrates a fun- 
damental link between the growth of supermassive BH and bulge formation. The connection 
between SFR and accretion rate has been explored by some authors (e.g., Heckman et al. 
2004; Hao et al. 2005; Satyapal et al. 2005; Dudik et al. 2005). Using the SDSS (Sloan Dig- 
ital Sky Survey) observations of 123,000 low-redshift galaxies, Heckman et al. (2004) found 
that the global volume-averaged SFR/M ratio is approximately 1000 in bulge-dominated 
systems, which is in agreement with the ratio of bulge to BH mass implied by the M bh — a 
relation (Marconi & Hunt 2003). Hao et al. (2005) also yielded similar results for several 
ten AGNs (including infrared-selected QSOs, optically-selected QSOs, NLSls). It was found 
that the ratio SFR/M of LINERs is apparently different from that of normal bright AGNs, 
of which the mass accretion rate M is estimated by assuming a constant radiative efficiency 
for all sources. The mass accretion rates M derived for LINERs may probably be incorrect 
(Satyapal et al. 2005), because RIAFs are believed to be in these sources and their radiative 
efficiencies are lower than these for standard thin disks. 

In this paper, we explore whether the LINERs follow the same correlation between 
accretion rate and SFR defined by normal bright AGNs. The accretion rates are re-estimated 
for these LINERs based on our spectral calculations for RIAFs. 
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2. The sample 

We adopt the sample given by Satyapal et al. (2005) for our present investigation, 
which includes 33 LINERs with detected hard X-ray nuclear emission (they are defined as 
AGN-LINERs in Satyapal et al. (2005)), 54 Seyferts, 15 quasars, 14 radio galaxies, and 
14 narrow-line Seyfert Is (NLSls). Several LINERs are radio-loud. All sources in this 
sample have estimated BH masses, FIR and bolometric luminosities. The BH masses of 
these sources are estimated through resolved stellar kinematics, reverberation mapping, or 
applying the correlation between optical bulge luminosity and central BH mass derived from 
nearby galaxies (see Satyapal et al. (2005), for more details). A significant fraction of galaxies 
in the sample are bulge-dominated and the SFRs of the host galaxies are calculated from 
Lfir by using the calibrated formula, in which the old stars' contribution has been properly 
subtracted (See Eqs. 5 and 6 in Satyapal et al. 2005). The calibrated formula adopted by 
Satyapal et al. (2005) can also be applied to some early type spiral galaxies in the sample, 
as this SFR-L FIR calibration is suggested to be applicable for all type galaxies (Kewley et 
al. 2002). All data can be found in Satyapal et al. (2005). We note that NGC 3215 is a 
Wolf-Rayet starburst galaxy, which is excluded in our present investigation. Some typos of 
their paper are fixed in this paper (see, Satyapal et al. 2006, for erratum). 



3. RIAF spectra 

In order to calculate the spectrum of a RIAF, we need to obtain the RIAF structure. We 
calculate the global structure of an accretion flow surrounding a Schwarzschild BH using a 
similar approach as Narayan, Kato, & Honma (1997). We adopt a simple a- viscosity (r r ^ = 
ap) and all radiative processes (Synchrotron, Bremsstrahlung and Compton scattering) are 
included consistently in our calculations for RIAF structure. The advection by ions and 
electrons has been considered in the energy equation, and a more realistic state of accreting 
gas (instead of a polytropic index 7 g ) is employed in the calculations, which is similar to 
that used by Manmoto (2000). The RIAF structure can be calculated when suitable outer 
boundary conditions are supplied. After the RIAF structure is obtained, the spectrum of 
the flow can be calculated in the same way as some previous authors (see, Manmoto 2000; 
Esin et al. 1996). 
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4. 



Results 



4.1. 



The relation between L 



X, 2-10 keV 



and rh for RIAF 



Three-dimensional MHD simulations suggest that the viscosity parameter a in the ac- 
cretion flow is ~ 0.05 — 0.2 (Hawley & Balbus 2002). Recent observations of black hole 
binaries (Maccarone 2003) and radio- loud AGNs (Marchesini, Celotti, & Ferrarese 2004), 
as well as theoretical accretion models (e.g., Narayan & Yi 1995; Esin, McClintock, & 
Narayan 1997) suggest that the critical accretion rate rh cr[t ~ 0.01, below which a RIAF 
is present. Our numerical calculations show that rh crit ~ 0.01 requires a ~ 0.2, and we 
adopt a = 0.2 in all our calculations. The parameter (3 (ratio of gas to total pressure in the 
accretion flow, (3 = P g /P tot ) is not an independent parameter, which can be related with a 
as (3 (6a — 3) /(4a — 3), suggested by MHD simulations (Hawley, Gammie, & Balbus 1996; 
Narayan et al. 1998). We adopt (3 ~ 0.8 for a = 0.2. Another parameter in our calculations 
is the outer radius i? ou t of the RIAF. We adopt i? ou t = 200i?s in all our calculations, where 
the Schwarzschild radius i?s = 2GM\^Jc 2 and is the BH mass. How the outer radius 
R out evolves with rh is still unclear, though it may be a function of rh. We find that the hard 
X-ray emission is mainly emitted from the inner region of the flow close to the black hole, 
which is insensitive to the value of R out adopted. A typical value 5 = 0.1 (the fraction of the 
dissipated energy directly heating the electrons) is adopted in all our calculations. Now, we 
can calculate the structure of the RIAF, and then its spectrum varying with accretion rate 
M when M hh is specified. 

For a fixed BH mass, = 1O 8 M , our calculation on the relation between 2-10 keV 
X-ray luminosity Lx, 2-10 kcv and accretion rate (from rh = 10~ 5 up to rh ~ m C rit ~ 10~ 2 ) 
is plotted in Figure 1. The dependence of 2-10 keV luminosity on rh can be roughly fitted 
by single power-law (solid line in Fig. 1), 



For rh < 10~ 3 ' 75 , Lx, 2-10 kcv °c ^ 2 and the 2-10 keV X-ray emission is mainly dom- 
inated by bremsstrahlung radiation. At higher accretion rate, 10~ 3 - 75 < rh < 10~ 325 , 
Lx, 2-10 kcv ^ 3 ' 5 an d both the Compton and bremsstrahlung radiation is contributes 
to the 2-10 keV X-ray emission. The Comptonized component will dominate the X-ray 
emission at higher accretion rate rh > 10~ 3 ' 25 which have L x , 2-10 kcv oc rh 2 ' 2 (see dashed 
line in Fig. 1). 

Our spectral calculations show that the 2-10 keV X-ray luminosity is roughly propor- 
tional to the BH mass. So we can calculate the accretion rate using the equation (1), if the 
BH mass and L x , 2-10 keV is given. 



\ogL x , 2-10 kcv = 2.37(±0.06) logm + 46.33(±0.22), 



(1) 
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4.2. The mass accretion rates of LINERs 

The central BH masses for all sources in this sample are available, so we can re-estimate 
the accretion rates for these 30 low- luminosity LINERs from their hard X-ray luminosity 
Lx, 2-10 kev based on our RIAF spectral calculations. The remaining three LINERs having 
X-ray luminosities higher than the maximal values expected by the RIAF model calculation. 
We estimate the accretion rates of these three sources through m = Lboi/^Edd, where Lboi 
is derived from the X-ray luminosity by using Lboi = 34Lx, 2-10 kcv- The Eddington ratio 
of 8 sources in the 14 radio galaxies are lower than 10~ 2 and these sources may be also in 
RIAF state. Four of these 8 sources (3C 31, 3C 84, 3C 338, 3C 465) have intrinsic X-ray 
luminosity (Evans et al. 2005) and the X-ray emission data for the remaining four sources 
(3C 88, 3C 285, 3C 327, 3C 402) are not available in the literature. We convert optical core 
emission measured by Hubble Space Telescope (HST) to a corresponding X-ray luminosity 
assuming spectral index a ox = 0.6 for these four sources (Marchesini, Celotti, & Ferrarese 
2004; Falcke, Kording, & Markoff 2004). Following the same way as LINERs, we re-estimate 
the accretion rates of these radio galaxies. The re-calculated mass accretion rates and the 
Satyapal et al. (2005) 's results are listed in Table 1 for comparison. 

In Satyapal et al. (2005), they showed that the distributions of SFR/M for LINERs and 
radio galaxies are apparently different from that of normal bright AGNs. Their estimates on 
M were carried out by assuming a constant radiative efficiency, which may be incorrect for 
low-luminosity LINERs with RIAFs. We plot the distributions of the ratio of SFR to the 
re-estimated mass accretion rate M for LINERs and radio galaxies in Figure 2 (upper and 
middle panels). We find that the distributions of SFR/M for LINERs and radio galaxies are 
similar to that of the normal bright AGNs in this sample (Seyfert, Quasar and NLS1, solid 
lines in Fig. 2). In Figure 3, we plot the relation of SFR and M for the whole sample using 
our estimates on mass accretion rates M. The linear regression yields for the normal bright 
AGNs (Seyfert, Quasar and NLS1): 

log SFR = (0.79 ±0.04) log M+ (1.66 ±0.06), (2) 

where SFR and M are both in units of M yr _1 . The linear regression for the whole sample 
is: 

\og SFR = (0.73 ±0.04) log M+ (1.55 ±0.07). (3) 

We find that the correlations are similar either for the whole sample or the subsample of 
normal bright AGNs (see, Fig. 3). We have not found any radio galaxies in this sample 
deviating significantly from the correlation. 
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5. Discussion 

The bolometric luminosities of AGNs can be estimated by integrating their spectral 
energy distribution or using an empirical relation (e.g., Woo, & Urry 2002; Elvis et al. 1994; 
Ho 1999), and then their mass accretion rates are derived if the radiative efficiency is known. 
A linear correlation between SFR and M was found by Satyapal et al. (2005) for normal 
bright AGNs, while the LINERs in their sample obviously deviate from this correlation, if 
the same constant radiative efficiency 0.1 is adopted for the whole sample. It was found that 
the accretion rates for LINERs are significantly lower than those of normal bright AGNs 
for given SFRs (see Fig. 10 in Satyapal et al. 2005). RIAFs are believed present in these 
low-luminosity LINERs, and their radiative efficiency e r is usually much lower than that 
of standard disks (Narayan & Yi 1994). This implies the accretion rates of LINERs are 
underestimated by Satyapal et al. (2005). In this paper, we derive the mass accretion rates 
of these LINERs from their X-ray luminosities based on RIAF spectral calculations, which 
should be more reliable. 

The standard thin disk is believed present in Seyferts and QSOs and its radiative ef- 
ficiency ~ 0.1. The slim disk may be present in the NLSls and its radiative efficiency 
should be also lower than that of standard disk due to the photon-trapping effects (e.g., 
Abramowicz et al. 1988; Ohsuga et al. 2002; Chen & Wang 2004). However, all the Edding- 
ton ratios of NLSls in our sample are less than 6 (£ Bol /L Edd , see table 2 in Satyapal et al. 
2005) and the radiative efficiencies of slim disks with such Eddington ratios may not deviate 
much from that of the standard disk (see, dashed line and solid line of Fig. 1 in Ohsuga et 
al. 2002). So, the derived mass accretion rates should be accurate to a factor of 2 for NLSls 
in these sample, which will not alter our main conclusion. 

The SFRs is calculated from the FIR luminosity for the whole sample since the FIR 
luminosity is widely used as a tracer of SFR in galaxies (e.g., Kennicutt 1998; Kewley et al. 
2002; Satyapal et al. 2005), though we can not rule out some contribution of dust heating by 
the AGNs and old stars in elliptical galaxies. However, there are lack of correlations between 
the Mid-IR and FIR for well-studied PG quasars, which suggests that they are dominated by 
different heating sources(e.g., Haas et al. 1999). If AGN heating dominates the cooler FIR- 
emitting dust, then there should be a correlation between quasar OUV and FIR luminosity, 
while such a correlation has not been found(e.g., McMahon, et al. 1999; Isaak et al. 2002; 
Priddey et al. 2003). In addition, QSOs and Seyferts roughly follow the same universal 
correlation between the FIR and radio emission deduced from "normal" galaxies, which 
strongly suggests that FIR emission is still powered by star formation rather by AGNs (e.g., 
Norris et al. 1988; Luo & Wu 2005). Although the FIR emission contributed by old stars 
in some elliptical galaxies could be important (e.g., very few sources in Bell (2003) sample 
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have old star contribution, about 4 of 249, higher 80 percent, see Fig. 6 in that paper), the 
average contribution of the old stars is around 30 percent. In our present SFR calculations, 
we use the calibrated formula, in which the old star contribution has been subtracted (Bell 
2003). So, in statistical sense, the SFRs estimated in our paper are reliable, though we 
cannot rule out the SFRs of a very few sources in this sample have been overestimated 3-4 
times. Even in this case, we believe this has not affected our main statistical results. 

We plot our re-estimated accretion rate M in Figure 3, and find that all LINERs follow 
the same correlation defined by the normal bright AGNs (e.g., Seyfert, Quasar and NLS1). 
The linear correlation between SFR and M indicates that black hole accretion evolves in the 
same way as star formation, which are both regulated by the interstellar gas in the host galax- 
ies. Our results show that the gases feed the black hole and star formation in low luminosity 
LINERs follow the similar way as luminous normal AGNs. Recent numerical simulations on 
the quasar activity triggered by the galaxy merger indeed show that the accretion rate M and 
SFR decreases simultaneously with time over several orders of magnitude after the quasar 
shines at its Eddington limit (Springel et al. 2004; Di Matteo, Springel, & Hernquist 2005). 
In their simulations, the gas in the galaxies is driven away by the bright quasar radiation, 
and then both the star formation and BH accretion are quenched simultaneously. The SFR is 
nearly linearly varying with M in nearly three orders (e.g., M ~ 5 x 10~ 4 — 5 x 10 _1 M yr _1 , 
SFR~ 5 x 10~ 3 — 5 Moyr^ 1 corresponding to ~ 1.7 — 2.0 Gyr) for the model with galaxies of 
virial velocity V vir = 160km/s (middle red lines in Fig. 2 of Di Matteo, Springel, & Hernquist 
2005). Our results are roughly consistent with their simulations, though the accretion mode 
transition has not been considered in their simulations. 

There are 8 radio galaxies (Eddington ratio less than 10" 2 ) in 14 radio galaxies, of which 
the accretion rates are derived from their X-ray emission. It is still debating whether the 
X-ray emission are dominated by the jet emission for some radio galaxies (Falcke, Kording, 
& Markoff 2004; Yuan & Cui 2005). If the jet emission in X-ray bands is important, the 
accretion rates derived in this paper are only the upper limits, and those 8 triangles in 
Figure 3 should be shifted towards left direction. As discussed in Sect. 4.2, we have not 
found these 8 sources deviating from the correlation between SFR and M, which may imply 
X-ray emission from the jets is unimportant at least for these 8 sources. 

We find that our results are insensitive to any model parameters except 5, which de- 
scribes the fraction of the viscously dissipated energy directly heating the electrons in the 
accretion flow. In our calculations, 5 — 0.1 is adopted, which is a typical value successfully 
used to model the observed spectra for some low-luminosity AGNs (e.g., Quataert & Narayan 
1999). Although winds may be present in RIAFs, the detailed physics is still unclear and 
such winds are only described by an artificial power-law parameter p w (e.g., Blandford & 
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Begelman 1999). In RIAF spectral calculations for X-ray wavebands, the value 5 is some- 
what degenerate with p w (e.g., Quataert & Narayan 1999). Our calculations show that 
similar conclusion can be obtained if S = 0.3 and p w = 0.9 are adopted for RIAF with winds. 

The present sample is a mixture of AGNs, of which both BH masses and bolometric 
luminosities are estimated. It is difficult to evaluate to what extent the sample is affected 
by the selection effects. A more robust sample is desired for testing this relation between M 
and SFR. 

The authors are grateful to the referee for constructive suggestions on our paper. We 
would like to thank Jufu Lu, Feng Yuan, Weimin Gu, Shuangliang Li and Yiqing Lin for 
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Table 1. Re-estimated accretion rates of LINERs and Radio galaxies 



Galaxy Name log 10 M a log 10 M b log 10 SFR c 

LINERs d 



NGC 4350 


-5.35 


-2.94 


-0.83 


NGC 1052 


-2.55 


-1.39 


-0.24 


NGC 3031 


-4.05 


-2.33 


-0.76 


NGC 4278 


-4.15 


-1.73 


-0.69 


NGC 4486 


-3.65 


-1.37 


-0.83 


NGC 4579 


-3.25 


-2.04 


0.09 


NGC 6500 


-3.95 


-1.89 


0.01 


NGC 4203 


-4.15 


-2.31 


-0.69 


NGC 4494 


-5.35 


-3.15 


-1.22 


NGC 4594 


-4.05 


-1.80 


-0.32 


NGC 4527 


-5.45 


-2.77 


0.44 


NGC 4125 


-5.55 


-2.61 


-0.62 


NGC 4374 


-5.15 


-1.99 


-0.69 


NGC 4696 


-4.05 


-2.06 


-0.47 


CGCG 162-010 


-2.35 


-0.92 


1.09 


IC 1459 


-3.85 


-1.74 


-0.62 


NGC 2787 


-5.95 


-3.39 


-1.16 


NGC 315 


-3.00 


-1.14 


0.01 


NGC 2681 


-5.30 


-3.07 


-0.47 


NGC 3169 


-4.13 


-2.32 


0.35 


NGC 3245 


-5.06 


-2.56 


-0.16 


NGC 3718 


-3.78 


-2.15 


-0.54 


NGC 4258 


-4.60 


-2.94 


0.01 


NGC 4261 


-3.71 


-1.73 


-0.39 


NGC 4410A 


-3.08 


-1.38 


0.27 


NGC 4457 


-5.00 


-3.21 


0.01 


NGC 4552 


-4.67 


-2.33 


-1.03 


NGC 4565 


-5.13 


-3.12 


-0.16 


NGC 6482 


-3.75 


-1.75 


-0.47 


3C 218 


-2.40 


-1.07 


0.72 


3C Radio galaxy d 








3C 31 


-3.751 


-2.12 


0.35 


3C 84 


-1.751 


-0.53 


1.18 


3C 338 


-3.451 


-1.64 


0.35 


3C 465 


-3.151 


-1.21 


0.35 


3C 88 


-3.351 


-1.41 


0.44 


3C 285 


-4.251 


-1.79 


1.15 


3C 327 


-3.351 


-1.10 


1.54 


3C 402 


-3.351 


-1.53 


0.44 



a Rcsult of Satyapal ct al. (2005), in unit of MQyr" 1 
b Our result, in unit of Moyr -1 
c in unit of MQyr -1 
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d The accretion rates of remaining 3 LINERs and 6 radio galax- 
ies are estimated from fa = ibol/^Edd with radiative efficiency 
ej! = 0.1. Because the Lx 2-10 koV of these sources are larger 
than the maximal values predicted by RIAF model with param- 
eters adopted in our paper and RIAFs is assumed not present in 
these sources. 
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-5 -4 -3 -2 

log m 



Fig. 1. — The filled circles connected by the dashed line show the integrated 2-10 keV X-ray 
luminosity versus the dimensionless accretion rate m predicted by the RIAF model. The 
solid line is a linear fit of L x , 2-10 keV to the model calculated points (filled circles). Model 
calculations are carried out for a 1O 8 M black hole, with a = 0.2, (3 = 0.8 and 5 = 0.1. 
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Fig. 2. — The histogram of SFR/M for the sources in the sample: 33 LINERs (the upper 
panel), 14 radio galaxies (the middle panel), and 81 normal bright AGNs (52 Seyferts, 15 
quasars and 14 NLSls, the lower panel). The dashed lines in upper and middle panels are 
results of Satyapal et al. (2005), while solid lines are ours. 
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Fig. 3. — Relation between SFR and M for the whole sample. The green-open squares are 
results of Satyapal et al. (2005), while the red-solid squares are for the re-estimated accretion 
rates. The red-solid line represents the linear fit to the whole 128 sources, while black-dashed 
line is for the normal bright AGNs. 



